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CNj [ Abstract 

A method allowing studies of the hadronic matter at the early evolution stage in A+A 

1 ^ ■ collisions is developed. It is based on an interferometry analysis of approximately conserved 
values such as the averaged phase-space density (APSD) and the specific entropy of thermal 
pions. The plateau found in the APSD behavior vs collision energy at SPS is associated, 
apparently, with the deconfincmcnt phase transition at low SPS energies; a saturation of this 

quantity at the RHIC energies indicates the limiting Hagedorn temperature for hadronic matter. 
It is shown that if the cubic power of effective temperature of pion transverse spectra grows with 
energy similarly to the rapidity density (that is roughly consistent with experimental data), then 
the interferometry volume is inverse proportional to the pion APSD that is about a constant 
ly-j ■ because of limiting Hagedorn temperature. This sheds light on the HBT puzzle. 

o: , , 

l/"") , Bogolyubov Institute for Theoretical Physics, Kiev 0314-3, Metrologichna 14b, Ukraine. 

PACS: 24. 10. Pa, 25.75.-q, 25.75. Gz, 25.75.Nq 
Keywords: relativistic heavy ion collisions, phase-space density, HBT correlations. 



o 

(N 



O 
=5 



X 



1 Introduction 



The main goal of experiments with ultra-relativistic heavy ion collisions is to study the new forms of 
strongly interacting matter which can be created under the extreme conditions. High densities and 
temperatures that arise in quasi-macroscopic systems formed during collision processes can result 
in the phase transitions from hadronic gas (HG) to Quark-Gluon Plasma (QGP) or sQGP pQ; the 
initial very dense pre-thermal stage of the collisions is, apparently, associated at RHIC energies with 
a specific form of matter - Color Glass Condensate (CGC) [2]. The bulk of hadronic observables 
("soft physics") are related, however, only to the very last period of the matter evolution, so 
called thermal or kinetic freeze-out - the end of the collective expansion of hadronic gas when the 
system decays. The evolution of wave function of a single particle, e.g. pion which can be created, 
annihilated and scattered many times in a dense stochastic (incoherent) surrounding, can be hardly 
considered. In fact, the single particle spectra bring information about the state of matter, e.g., its 
temperature and flows, at the very end of the hadron gas evolution, and the HBT interferometry [Sj 
is related directly to the structure of emission function, i.e. the space-time density of last hadronic 
collisions. Thus, the correlation measurements itself cannot be related directly to the preceding 
hot and dense stages of the matter evolution in A+A collisions where a formation of new forms 
of matter is expected. Our basic idea is to use the "conserved observables" which are specific 
functionals of spectra and correlation functions - integrals of motion. Then such an observables, 
being conserved during the matter evolution (or some period of the evolution) can be related to 
the state of matter at the early evolution stage. 
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A structure of the integrals of motion, besides the trivial ones associated with the energy, 
momentum and charges, depends strongly on a scenario of the matter evolution. The actual 
numbers of partonic degrees of freedom released in ultra-relativistic nuclear-nucleus collisions are 
fairly big, up to tens of thousands. As it is commonly supposed such a quasi-macroscopic system 
becomes thermal at the early stage of the collision process: the time of thermalization for the 
RHIC energies is estimated to be from 0.6 fm/c (the phenomenology analysis [I]) to 1 fm/c (the 
"black hole thermalization" [2]) and even more, 3 fm/c (the pure pQCD result [2]). Then the system 
expands nearly isoentropically - the latter is standard point in the hydrodynamic approach to A+A 
collisions and is advocated theoretically for the RHIC energies ^Ej. As for the SPS energies, the 
approximation of perfect hydrodynamics successfully describes the spectra and correlations while 
overestimates the elliptic flows [Hj. The latter can be connected with some viscosity effects at this 
energies. Keeping this in mind, one can, nevertheless, consider, at least approximately, the entropy 
produced in A+A collisions as integral of motion that carries out the information as for very initial 
thermal stage of these processes. The another "conserved observable" is found recently it 
is the pion phase-space density (/) averaged over momentum (totally or at fixed rapidity) and 
configuration space. It is also about a constant during the stage of chemically frozen expansion. 
So, measurements of the entropy and average phase-space density (APSD) in thermal hadronic 
systems at the final, freeze-out stage of A+A collisions makes it possible to look into the previous 
stages such as the partonic thermalization and the hadronization (or chemical freeze-out). 

In the paper we express the entropy and APSD of thermal pions through the observed spectra 
and interferometry radii irrespective of unknown form of freeze-out (isothermal) hypersurface and 
transverse flows developed based on the approach proposed in the Ref. [OJ- Our aim is to study the 
properties of the matter around hadronization stage at different energies of A+A collisions, from 
AGS to RHIC, provide the general analysis of the results and to make conclusions as for possible 
new forms of matter formed in these processes. 

2 The entropy and APSD as observables in A+A collisions 

To clarify the problem let us start from the entropy of thermal pions. As well known the expression 
for entropy of a gas of bosons/fermions has the following covariant form 

S=(2J + 1) J ^0 [-(2^) 3 /ln((2vr) 3 /) ± (1 ± (2vr) 3 /)m(l ± (2^) 3 /)], (1) 

where ± sign corresponds to bosons/fermions with total spin J and a is some hypersurface in 
Minkowski space. The value depends on the distribution function, or the phase-space density 
f(x, p), that should be known to make the corresponding estimates. It is easy to show, however, that 
the phase-space density, e.g. of tt~ , cannot be extracted from the two- and many- particle spectra 
even if the system at kinetic freeze-out is characterized by the locally equilibrated distribution 
function. To make it clear let us write the Wigner function (an analogy of the phase-space density 
for the quantum systems) for weekly interacting particles in the mass-shell approximation |12j : 

f(x,p) = (2vr)- 3 J d i q5(q-p)e-^ (a + (p - (l/2)g)a(p + (1/%)) . (2) 

1 Unlike many authors (see, e.g., |1UI Ml] ) we utilize almost totally averaged value of the APSD, {/), instead of 
momentum-dependent one (f(p)) because the extracted value of the later quantity is affected essentially by assumed 
form of freeze-out hypersurface as well as a assumed profile and intensity of flow on it, see Ref. and correspondent 
discussion in Ref. 
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Here the brackets (...) mean the averaging of the product of creation and annihilation operators 
with a density matrix referred to the space-like hypersurfaces where particles become or are already 
nearly free. If the freeze-out is sudden, one uses usually a thermal density matrix at the freeze-out 
hypersurface. The invariant single- and double- (identical) particle spectra have the forms: 

n (p) = (a + (p)a(p)) ,n(p 1 ,p 2 ) = n(p 1 )n(p 2 ) + \(a + \pi)a{jp 2 ))\ 2 , (3) 

and the correlation function is defined as the following 

C(p, q) = n(p 1 ,p 2 )/n(p 1 )n(p 2 ). (4) 

It is easy to see from Eq. © that the complex phase of the two- operator average (a + (pi)a(p 2 )) 
cannot be extracted from the single- and double- particle spectra. It is possible to show that the 
same takes place even when many-particle spectra are included into an analysis. Therefore one 
cannot reconstruct the distribution function f(x,p), see Eq. (J2J), in a model independent way. 
However, since f{x,p) is real (note, that for the locally equilibrated quantum systems the Wigner 
function is positive), one can use Eq. © to express the f 2 (x,p) integrated over space coordinates 
just through squared absolute value of the two- operator average, \(a + (p — (l/2)q)a(p + (l/2)q))f, 
integrated over q. Then, accounting for the direct links (j2J), @ of the latter value with the 
single particle spectrum (that is just the f(x,p) integrated over space coordinates) and correlation 
function, one can get the phase-space density averaged over some hypersurface a, where all particles 
are already free: a = a out, and over momentum at fixed particle rapidity, y = 0, 



_ f(f(x,p)) 2 p»da,d 2 PT _ (2^)- 3 J Po l n 2 (p)(C(p,q) - ljAjV 
V{<r,V))y~o f f(x,p)pHda^p T dN/dy 



directly from the experimental data in full accordance with the pioneer Bertsch idea JB]. Using 
the standard Bertsch-Pratt parametrization the last equality in (JSJ) can be re-written through the 
interferometry radii, as it is presented in Eq. 0. 

On the face of it, the extracted value of the APSD © does not help to determine the entropy l[T|). 
and so, some phenomenological functions that can reproduce the approximate Gaussian behavior 
of the correlation function are usually suggested ^3]. However, since we cannot extract the phase 
of the two- operator average, there is the infinite set of distribution functions compatible with the 
observables and, therefore, the entropy calculated will depend on the class of the functions we 
choose which. Here we propose the method to estimate the entropy using just the APSD without 
assuming any concrete expression for the phase-space density. 

The method is similar to what was proposed in Ref. [S] for analysis of the overpopulation of the 
phase-space. The idea is based on the standard approach for spectra formation JH] that supposes 
the thermal freeze-out in expanding (with 4-velocity field u u ) locally equilibrated system happens 
at some space-time hypersurface with uniform temperature T and particle number density n (or 
chemical potential //). Then, within this approximation which is probably appropriate in some 
"boost-invariant" mid-rapidity interval, Ay < 1, the phase-space density of pions, 

/ = fi, q .(x,p) = (2n)- 3 (eM Mx) ^~ fl ) ~ I)" 1 , (6) 

totally averaged over the hypersurface of thermal freeze-out, a = oth, and momentum except the 
longitudinal one (rapidity is fixed, e.g., y = 0) will be the same as the totally averaged phase-space 
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density in the static homogeneous Bose gas [Sj: 



fd 3 pl 2 2 ^ 5/2 f [roRsRl (2nm d T Zi T dy) ) dm T 

^(MrtU-i^- ^ ^ , P) 



where 



/^^(exp^o-M)-!) -1 (8) 

and (3 and coincide with the inverse of the temperature and chemical potential at the freeze-out 
hypersurface. 2 The k = 1 if one ignores the resonance decays. Here we neglect interferometry 
cross-terms since they are usually rather small in the mid-rapidity region for symmetric heavy ion 
central collisions at high energies. 

The last expression in Eq. (J7| just corresponds to the last term in © calculated in the Bertsch- 
Pratt parametrization. The important result presented by the first equality in is based on 
the properties of relativistic invariance of the distribution function © and its local isotropy in 
momentum in the rest frames of each fluid element. Then, using the "boost-invariance" within 
homogeneity length, Ay ~ 1, the integrals over p^da^cfpx of diverse functions Fr* of the locally- 
equilibrium distribution, F/^(fi_ eq _), contain the common factor, "effective volume" Veff — J ^ 
(r/ is rapidity of fluid), that completely absorbs the flows u^(x) and form of hypersurface cr(x) in 
mid-rapidity. For instance, if Fu\(fi. eq .) = fi. eq ., then 



^ = J p»da^p T f Leq . = V eff J d 3 p^ = n th V eff (9) 



where n t h is thermal density of equilibrium ideal Bose gas, and similar takes place for F^(fi eq ) = 
ff eq in ©• Thus, the effective volume F e // is cancelled in the corresponding ratios. This factor- 
ization property has been found first for Eq. Q in Refs. |16L 117) . multiple used for an analysis of 
particle number ratios (see, e.g., Ref. JH]) and recently generalized for a study of the APSD in 
Ref. 0. 

In this work we apply the approach to an analysis of the thermal pion entropy per particle, 
or specific entropy of thermal pions. Using the same approximation of the uniform freeze-out 
temperature and density and Eq. Q with some local equilibrium distribution we get the following 
expression for specific entropy in mid-rapidity: 



dS/dy = f d 3 p [-f eq In f eq + (1 + f eq ) ln(l + f eq ) 
dN/dy $d?p-f eq 



(10) 



In the above ratio due to the factorization property the effective volume is cancelled and the final 
expression depends only on the two parameters: the temperature and chemical potential at freeze- 
out. The temperature can be obtained from the fit of the transverse spectra for different particle 
species and we will use the value T = 120 MeV as a typical "average value" for SPS and RHIC 
experiments. 



2 Note here that (27r) 3 / is unitless not only in the natural units where H = c = 1 but also in conventional system 
of units (e.g. SI) where (27r) 3 / is expressed via the Plank constant: (2irh) 3 f = h 3 f. Therefore the above value can 
be interpreted as distribution of the population numbers in elementary phase-space cells unlike phase-space density 
/ that carries units h~ 3 . 
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Another parameter, the chemical potential, cannot be extracted from the spectra, its value 
could be fairly high even for the thermal pions because of chemical freeze-out and this parameter 
is crucial for an estimate of the entropy. We extract the chemical potential from an analysis of the 
APSD following to l[7|). The factor k is accounting for a contribution of the short-lived resonances 
to the spectra and interferometry radii and absorbs also the effect of suppression of the correlation 
function due to the long-lived resonances [S]. Because of the chemical freeze-out a big part of 
pions, about a half, are produced by the short-lived resonances after thermal freeze-out. It leads 
to an increase of the APSD despite the maximal particle and entropy densities of pions at post- 
hydrodynamic stage of the evolution is reached at the end of the hydrodynamic expansion - at 
thermal freeze-out as discussed in detail in Ref. [§]. To estimate the thermal characteristics and 
"conserved observables" at the final stage of hydrodynamic evolution by means of Eqs. IjTjl. (|10|) 
one needs to eliminate non-thermal contributions to the pion spectra and correlation functions from 
resonance decays at post freeze-out stage. To do this we use the results of Ref. [H| where a study 
of the corresponding contributions within hydrodynamic approach gives the values of parameter k 
to be k = 0.65 for SPS and k = 0.7 for RHIC, if half of pions is produced by the resonances at post 
freeze-out stage. Then, from Eqs. 0, (jHJ) one can extract the pion chemical potential at thermal 
freeze-out. This makes it possible to estimate the average phase-space density, the specific entropy 
of thermal pions and other thermal parameters of the system at the end of the hydrodynamic 
expansion as explained above. 

3 The analysis of experimental data and the results 

To evaluate the APSD of negative pions by means of Eq. we utilize the yields, transverse 
momentum spectra and interferometry radii of tt~ at mid-rapidity measured in central heavy ions 
collisions by the E895 and E802 Collaborations for AGS energies HSIEU, NA49 Collaboration 
for SPS CERN energies HH H21 H31, STAR and PHENIX Collaborations for RHIC BNL energies 

To calculate the APSD we have to integrate over the whole pt region. Always, if possible, we use 
the pion yields and values of radii in each measured pt bin instead of analytical parameterizations 
for the transverse spectra and interferometry radii. To interpolate between successive data points we 
use an polynomial functions, the degree of the polynomial curves is chosen to be 2 or 3. Outside the 
whole measured pt region we make analytical extrapolations with parameters that catch the main 
tendencies of the data points in the measured pt region. We utilize the analytical parameterizations 
that are typically used in fitting procedures and are motivated by hydro model calculations. Namely, 
for the transverse spectra at high pt we assume the exponential parameterizations with slopes that 
correspond to average ones in the measured px region (that is, actually, the average slope in the 
measured pt < 1 GeV points). For small unmeasured pt our exponential extrapolations are taken 
with the slopes which are defined from a requirement of coincidence of the resulting total 7r~ 
yields at mid-rapidity with the values presented by experimental Collaborations. To extrapolate 
interferometry radii behind the measured px bins, we utilize the widely used simple analytical 
parameterizations for interferometry radii, a/(b + c * niT) d , with numerical parameters which are 
taken, if possible, from the fit given by the experimental Collaboration, or we determine them 
ourselves if they are not presented. Also, we use the approximation Ro = Rs for the outward and 
sideward interferometry radii below the minimal measured pt momentum, because this equality 
should take place for px = due to evident (and well known) symmetry reasons. We neglect 
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interferometry cross-terms because they are rather small at mid-rapidity (see, e.g., |23j)). and, for 

example, at SPS energies yj ' R? R\ — Rq l is equal to RoRl with a few percents accuracy. Since 

the interferometry radii are measured by PHENIX Collaboration for — 30 % centrality events at 
y / s/viv = 200 GeV |26| , we increase the interferometry volume measured by PHENIX Collaboration 
at this cm. energy by a factor of 1.215 to get the interferometry volume corresponding to the most 
central — 5 % centrality bin in accordance with N part dependence of the Bertsch-Pratt radius 
parameters found in Ref. |26| . 

The results for the APSD at mid-rapidity for all negative pions at the AGS, SPS, RHIC energies 
in logarithmic cm. energy scale are presented in Fig. 1. Note that for SPS energy domain we 
use the results of NA49 Collaboration since it presents interferometry data as well as the data 
on transverse momentum spectra which are necessary to get the APSD values, while the CERES 
Collaboration demonstrates the interferometry radii only. 

The non-monotonic structure of the APSD behavior seen in Fig. 1 in the AGS-SPS energy 
domain could be an indicator of new physical phenomena as it is discussed in next Section, so a 
statistical significance of this structure is important and we demonstrate it in Fig. 2 at AGS, SPS 
energies in ^/snn scale. The comment is the following. Because of fast (exponential) decrease in of 
transverse momentum spectra with the uncertainties of the spectra and HBT radii in the region 
of high pt has only a little influence on the calculated APSD. The latter are mostly affected by the 
values of the HBT radii at low pj>. The most important are systematic (from a choice of analytic pr 
parameterizations) uncertainties in this region behind the lower measured pt bin, and experimental 
errors in the HBT data. Taking into account that the experimental errors in the HBT radii do not 
exceed typically 0.5 fm (see, e.g., |23|). and estimating the resulting systematic uncertainties from 
our choice of pt parameterizations as 0.5 fm, we get the uncertainty of 1 fm in each HBT radius. 
For a rough estimate of relative errors in the APSD values let us assume that average value of HBT 
radius in the region of interest is approximately 5 fm. Then, supposing that all uncertainties are 
independent and accidental, we can estimate that those errors are approximately 30 percents. As 
one can see from Fig. 2, the APSD in AGS energy domain can be approximated with quite good 
accuracy by linear in ^Jsnn function, cl^/snn + b, where a = 0.03768 GeV _1//2 and b = 0.01345. 
Then the extrapolation of this tendency to SPS energy domain results in rather high values for 
APSD, for example, at y / J/v r /v = 12.3 GeV we get (27r) 3 (/) = 0.477 that is far from the estimated 
error bars, (2-7r) 3 (/) = 0.203 ± 0.06. While the final conclusion about statistical significance of the 
observed tendencies needs in more detail studies both theoretical and experimental, our estimate 
from the above analysis is that it is very likely that experimental data really indicate non-monotonic 
behavior of the APSD as a function of collision energy in the AGS-SPS energy domain. 

In Fig. 3 the APSD of thermal negative pions at the SPS and RHIC energies are presented. 
Here and below we demonstrate the values for the thermal negative pions at the RHIC energies 
which are mean values of STAR and PHENIX data. Note, that in Fig. 3 we present also the APSD 
of negative thermal pions at chemical freeze-out in assumption of chemical equilibrium, (p, n = 0). 
We use the temperatures of chemical freeze-out for different colliding energies from Refs. |27j . 
where they were found from analysis of particle number ratios. We do not calculate the APSD of 
thermal pions at the AGS energies because estimates of resonance contributions to the pion spectra 
have been done in 9 for the SPS and RHIC energies only. 

The APSD of negative thermal pions are used then to extract the chemical potentials [i of them 
at thermal freeze-out with = 120 MeV at different SPS and RHIC energies, and after that 
to calculate the specific entropies, s = ^ju y lfTo"|) . the entropies, ^| = s * and the densities 
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nth, see Eq. ©, of negative thermal pions. The values founded are presented in Figs. 4, 6, 7. 
Also we demonstrate in Fig. 5 the interferometry volume Vi n t = (2vr) 3//2 i?o^5^L calculated at 
small pt — 0.06 — 0.07 GeV as the function of the rapidity densities, dN* /dy, of all negative 
pions at mid-rapidity in central nucleus-nucleus collisions. For the RHIC energies we use in Fig. 5 
the interferometry radii measured by STAR Collaboration because the Collaboration presents the 
interferometry measurements in lower pt bins as compare to PHENIX Collaboration. If there are 
no experimental data in selected pr bin, px — 0.06 — 0.07 GeV, we calculated the correspondent 
values using the analytical parameterizations of the interferometry radii as explained above. The 
obtained values, Vint-, are used then to evaluate the ratios (dN n /dy)/Vi n t that are demonstrated 
and compared with the thermal densities © in Fig. 6. 

In Fig. 8 we present, in addition to Fig. 7 where the entropies dS/dy of negative thermal pions 
are demonstrated, the rapidity densities, dN w /dy, of all negative pions at mid-rapidity in central 
nucleus- nucleus collisions for the AGS, SPS and RHIC energies. The experimental values for the 
pion rapidity densities are taken from Refs. ^§1 I2H1 Ell 1221 EH ESJ ■ We used the rapidity densities 
of pions instead of those for all charged particles, since pions are the new produced particles which 
are not contained initially in colliding nuclei, and, therefore, more directly represent a mechanism 
of the particle production in A+A collisions. It is especially important for collision processes 
with relatively low multiplicities at the AGS energies, where a large fraction of registered charged 
particles are protons which were not produced in collision processes being initially in colliding 
nuclei. The lines in Fig. 8 represent the logarithmic law of energy dependence for negative pion 
multiplicities: alog 10 (y / J/v r /v/6), where a = 160(230), b = 1.91 GeV (3 GeV) for solid (dashed) lines 
respectively. Note, that we use the STAR point dN™ /dy = 249 for ^snn = 130 GeV because this 
value - the result of a use of the Bose-Einstein distribution as a fit function, is closer then another 
STAR value based on a Boltzmann-like fit (see detail in [2l]) to the value dN w /dy = 270 ± 3.5 
reported by PHENIX Collaboration and is also closest to the result dN 77 /dy = 287 ± 20 that is 
deduced by STAR Collaboration from their measurements of negative hadrons, antiprotons and 
negative kaons spectra in Ref. |28j . Note that the latter value (it is not presented in Fig. 8) is even 
closer to the fitting line, which is showed in Fig. 8 as the solid line, than demonstrated experimental 
value. 

4 Discussion and interpretation of the results 

Let us start from an analysis of the ^/snn dependance of the averaged phase-space density per 
unit of rapidity (APSD). Figure 1 is based on Eq. © or the last term in Eq. with k = 1 
and demonstrates a behavior of the "raw" APSDs (which are, actually, the "asymptotic" APSDs 
related to the times when pions are detected) accounting for all negative pions, ir~ , thermal and 
from resonance decays at post thermal stage. One can see that the APSD grows significantly with 
energy at the AGS energies, then has the plateau starting from the lowest SPS energy, 20 AGeV, 
till 80 AGeV and then begins to grow again, apparently very slowly at RHIC as one can conclude 
from the non quite compatible experimental data of the STAR and PHENIX Collaborations. 

Unlike a fast decrease of particle n(x) and phase-space f(x,p) local densities, the totally aver- 
aged phase-space densities of thermal pions are conserved during a chemically frozen evolution [0] 
Roughly, it is proportional to the total ("raw") APSD at the SPS and RHIC energies, see Fig. 3. If 
the same properties take place at the AGS energies too, then one can easily interpret the behavior 
of the averaged phase-space density in Fig. 1. When the energy of collisions at AGS grows, the 
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initial hadronic density and phase-space density increase; since the pion APSD is conserved, its 
observed value also grows. Its rise stops at low SPS energies, this means that the initial density of 
pions also stops increasing. The simplest explanation is: an excess of the initial energy begins to 
be transformed into new, non pionic (hadronic) degrees of freedom, possible to quarks and gluons. 
The pure hadronic stage appears later when the densities became smaller than the initial ones, 
therefore the initial APSD of pions depends then not on the initial energy density but on the den- 
sity determined by the hadronization temperature T c . Figure 3 demonstrates that the APSD of 
the thermal pions grows indeed very slowly starting from high SPS energies that can reflect the 
fact of saturation of the temperature of the phase transition at the RHIC energies. The APSD at 
the thermal freeze-out is slightly higher than at chemical one. It is because the conservation of 
the APSD should take place at perfectly chemically frozen hadronic evolution; there is, however, a 
residual effect of increase of pion number because of an excess of resonance decays into the expand- 
ing gas over back processes of the recombination. This difference does not contradict to typical 
estimates that roughly 2/3 of pions at hadronization stage are " hidden" in resonances |29| 1301 118j 
and about half of pions are already thermal to the end of the hydro evolution. The pion APSDs 
at the chemical freeze-out are calculated using the thermal parameters of that stage that were 
found from an analysis of particle number ratios at SPS and RHIC in Refs. |27| . A dependance of 
the chemical potential fi of thermal pions at thermal freeze-out on y / ijv r /v is demonstrated in Fig. 
4. One can see that the chemical potential of thermal pions at freeze-out is saturated somewhere 
between 50 and 60 MeV. 

Figures 5 and 6 are related to the behavior of the interferometry volume Vi n t on multiplicities, 
dN w /dy, in central collisions at different energies of AGS, SPS, and RHIC. The interferometry 
volume has a tendency to grow over a broad energy range, as one can see from Fig. 5. However for 
central Pb+Pb (Au+Au) collisions at different energies, the corresponding increase is much slow 
than the proportionality law between Vi n t{sfs) and dN v /dy(^/s). Also, there is the statistically 
reliable violation of a monotonic behavior of V% n t |31j manifested in the evident decrease of Vi n t 
in AGS energy interval. The phenomenon has been associated with the supposed constancy of 
the (kinetic) freeze-out value of pion mean free path while the transition from the nucleon to pion 
dominated matter happens within that energy range Note that the relatively steep rise of 

the pion APSDs (see Figs. 1 and 2) at a moderate increase of pionic multiplicities at the AGS 
energies (see Fig. 8) is caused by the discussed decrease of Vi n t with beam energy. Of course, the 
conservation of the APSD in absence of the deconfinement phase transition leads in any case to 
some rise of its value with initial energy density. As follows from Fig. 6, the dN n /dy grows with 
energy significantly faster than Vint- This fact is the main component of the HBT puzzle |32l 133 j . 
To understand it qualitatively, let us very roughly estimate the APSD supposing that the 
transverse spectra have mainly exponential behavior vs transverse mass m-r, oc exp (— rriT/T e ff), 
where effective temperature T e ff depends on the thermodynamic temperature at the hypersurface 
of thermal freeze-out a and flows at a. Then, assuming that integral I over dimensionless variable 
m,T/T e ff depends on energies of collisions fairly smoothly, one can write 

vu-fi><*ig&- (ii) 

where the interferometry volume is taken here at the smallest tut- Thus, a proportionality between 
Vint and the particle numbers dN/dy is destroyed by a factor (f)T^jj. So, if the APSD and Vint 
only slightly grow with energy, mostly an increase of Tf// com d compensate a growth of dN/dy in 
Eq. (jllj). One can see that it is the case: for example, the ratio of cube of effective temperatures of 
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negative pions at y / i]v7V = 200 GeV (RHIC) to one at 40 AGeV (CERN SPS) gives approximately 
2, while the ratio of corresponding mid-rapidity densities is approximately equal to 3 |211 1241 05] . 
It could be only in the case of an increase of the pion flows in A+A collisions with energy. If 
the intensity of flows grows, it leads to a reduction of homogeneity lengths and the corresponding 
interferometry radii [31] • This effect almost "compensates" a growth of final geometrical sizes of 
the system with energy in observed interferometry volumes. Indeed, as one can see from Fig. 6 the 
freeze-out densities for pions, nth, become noticeably smaller than formally defined HBT densities, 
dy , starting from the high SPS energies. In other words, the interferometry volume at those 

Vint 

energies becomes to be significantly smaller than the effective one occupied by the system and 
defined by Eq. @. 

The result (fTTj) brings about some more details. If at any fixed energy y / s]v r /v we look at 
the evolution in time of Vint, we found that it can be nearly constant since the values dN/dy, 
APSD (/) and effective temperature T e ff in r.h.s. of Eq. are approximately conserved for 

the thermal pions during the chemically frozen hydro-evolution As the result, the "HBT 
microscope" at diverse energies "measures" the radii that are similar to the sizes of colliding nuclei. 
It provides an explanation to the phenomenological observations 23 1 3S] that in central Pb+Pb and 
Au+Au collisions the interferometry volumes grow rather slowly with energy, and only due to the 
longitudinal interferometry radius grows (transverse sizes are equal), while, at the same energy, the 
Vi n t depends strongly on the sizes of colliding nuclei and on the impact parameters in non-central 
collisions. 

Our other observations are based on one more conserved value, the entropy. In principle, the 
entropy S of the pion component alone can be changed even in the perfect fluid as well as the pion 
numbers N (see discussion above). One can expect, however, that such deviations will be small 
for the ratios s = dN/'dy or ^ or ^ e s P ec ifi c entropy. Then using the result (|1U|) for the latter value 
at the thermal freeze-out and the chemical potential extracted from an APSD analysis, one can 
determine the specific entropy of pions. The corresponding estimates give the values of s to be 
approximately equal to 4 at SPS except for the top SPS energy, where s ~ 3.79. At RHIC energies 
for ^snn = 130 GeV, the corresponding averaged value is equal to 3.66 (one can conclude from 
Fig. 1 that such a low value is, probably, artefact and result of relatively high discrepancy between 
the STAR and PHENIX data), and for ^snn = 200 GeV the specific entropy of negative pions is 
equal to 3.82. 

The total entropy of thermal negative pions per unit of rapidity, that is s multiplied by the 
rapidity densities of thermal pions, (dN w /dy)/2, is presented in Fig. 7. One can see that the 
entropy starts to grow faster at the SPS energies but at the RHIC energies this tendency is canceled. 
To understand better the situation let us look at the "raw" data in Fig. 8 representing a behavior 
of the negative pion rapidity density in the AGS, SPS and RHIC experiments. From this picture, 
that makes the tendencies which are seen in Fig. 7 more evident, one can conclude that at the SPS 
energies there is, indeed, an anomalously large slope of an increase of the pion entropy (and the 
number of pions) with energy. The observed multiplicities of negative pions at SPS at 80 AGeV 
(158 AGeV) are by factors of 1.09 (1.143) larger than ones extrapolated in accordance with the 
tendency (dashed line) observed at the AGS energies and, apparently, at the RHIC energies where 
relatively low pion multiplicities could be remnants of CGC formation 2 . Taking into account that 
chemically equilibrated thermal model correctly describes particle number ratios, including pions, 
at the AGS and RHIC energies [301 and that at the top SPS energy there is a problem of a "pion 
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deficit" we can suppose that the observed "excess" of pions is caused by a mechanism shifting 
the pion production at the SPS energies from the equilibrium. This mechanism could be also in 
some degree responsible for the reduction of K/tt ratios as compare to ones in chemical equilibrium 
model |37j - the effect was observed by NA49 Collaboration ("horn" puzzle) [22]. Some increase 
of the pion APSD at SPS 158 AGeV, see Figs. 1, 2, 3 could be also caused by the "extra pion" 
production at high SPS energies. 

What could be the reason of the "extra pion" production out of equilibrium at SPS energies? 
The intriguing possibility is that such an effect is the manifestation of the QCD critical end point 
(CEP) which is the terminating point of the first order phase transition line (about physics of 
the CEP and its location at QCD phase diagram see, e.g., [SHI)- Indeed, since the CEP acts 
as an attr actor of the isentropic trajectories of the thermodynamic system evolution [HH], the 
critical domain can influence the particle spectra for some range of collision energies, e.g., this 
could be responsible for the "step" behavior of the kaon effective temperature, discovered by NA49 
Collaboration at the SPS energies [22]- Then nonequilibrium features, which accompany the phase 
transition in the expanding systems, e.g., a rise of the bulk viscosity in the mixed phase due to 
a finite relaxation time and variation of sound velocity in the transition region [301 ' could lead to 
a dissipation of kinetic energy and to entropy production. Another effect, that can result in the 
"extra pion" production in A+A collisions near the CEP, is a significant reduction (that is maximal 
in the vicinity of the CEP at the QCD phase diagram) of a meson mass from its vacuum value, and 
mass shift of other resonances (p, etc.) due to sigma exchanges ^2E2]. As a result these species 
are rather numerous around the CEP. As m CT < 2m n in a vicinity of the CEP [411 142j , such sigma 
mesons cannot decay into nn in this state and can do it only after a meson masses are increased 
that happens when the density in A+A collisions are reduced by the system expansion. One can 
speculate that when it happens the rate of inelastic collisions can be not high enough to push 
the new produced pions into chemical equilibrium, while the elastic collisions can still thermalize 
them. That is why there could be no peculiarities in the pion transverse momentum spectra in low 
Pt region, as predicted in Ref. [l^, while an enhancement of pion yields can be considered as a 
possible manifestation of the CEP in A+A collisions. 4 Since the "pion excess" is maximal at the 
highest SPS energy 158 AGeV, it could mean, unlike present expectations , that the CEP is 
situated in the QCD phase diagram closer to the chemical freeze-out point at highest SPS energy 
than at the lowest SPS one. This possibility is argued for recently in Refs. [Sj and [3H]- It is 
noteworthy that because of an inhomogeneity of the baryonic chemical potential and temperature 
in rapidity at the chemical freeze-out hypersurface (see, e.g., an analysis that has been done in Ref. 
46 ), the condition for the thermodynamic trajectory of system evolution to pass around the CEP 
could be realized just in mid-rapidity region. Then the enhancement of pion production can be 
observed from diverse particle number ratios in the unit central rapidity interval rather than from 
47T abundances. 

At the RHIC energies both mechanisms, bulk viscosity and sigma mass reduction resulting in 
intensive entropy and pion multiplicities rise, can be inefficient since at small net baryon density, 
which are typical for that energies, the crossover far from critical domain around CEP might, 
apparently, happen (see, e.g., |17jb 

A few remarks are in order here. Presently, one of the interpretation of a larger rate of an 

3 The problem of pion deficit in mixed (full 4-7T geometry - midrapidity) particle number ratios in central 158 AGeV 
Pb+Pb collisions was pointed out earlier in Ref. |.'-Sfi| were preliminary data on hadron multiplicities were analyzed 
within hadron gas models. 

4 The "extra pion" production was considered earlier as the signature of the chiral phase transition in Ref. |43|. 
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increase of pion production in A+A collisions at the SPS energies as compared to AGS energies 
and properly normalized p + p(p) collisions is based on the statistical model of the early stage 
(SMES) ^43] (see also [12] )■ In this model the kink-like change (and horn-like structure of K + /ir + 
ratios) is a direct consequence of an onset of the deconfinement and liberation of massless quark- 
gluon degrees of freedom in the Lorentz-contr acted fireball at the initial stage of A+A collisions. 
However it seems that both basic assumptions, namely, massless quark-gluon plasma and the 
Lorentz-contracted fireball, are not supported by further studies. The effective masses of quarks 
and gluons in the QGP are temperature-dependent and grow with temperature (see, e.g., Ref. 
|50| where the validity of quasiparticle picture of the quark-gluon plasma at high temperatures is 
advocated) that turns down the possibility to treat the QGP as weak coupling massless quark- 
gluon system even for T > 3T C . It is also clearly seen from the lattice QCD results |51| where 
the pressure and energy density are both below the Stefan-Boltzmann limit even at very high 
temperatures which will be hardly reached in heavy ion collisions. Another crucial assumption in 
the SMES is a formation of the longitudinally Lorentz-contacted fireball in rest, that is actually the 
Landau- type initial conditions for hydrodynamic expansion. It was demonstrated [52], however, 
that the Landau-type initial conditions are unable to reproduce effective temperatures together with 
other data (multiplicities and rapidity distributions) at the SPS energies, and that these quantities 
can be described altogether only when one uses large initial volume with an appropriate velocity 
distribution (see also [53]). 

5 Conclusions 

The method allowing analysis of an early stage of the hadronic matter evolution in A+A collisions is 
developed. It is based on studies of approximate integrals of motions for the evolution of hadronic 
systems, such as the totally averaged phase-space density (APSD) and the specific entropy of 
thermal pions. We express these quantities through experimental data on the spectra and Bose- 
Einstein correlations in a way that does not depend on the freeze-out hypersurface and collective 
flows developed. Our estimates of the APSD at hadronization stage are close to the corresponding 
ones that we calculate based on the results of analysis of particle number ratios. A behavior of 
the pion APSD vs collision energy has a plateau at low SPS energies that indicates, apparently, 
the transformation of initial energy to non-hadronic forms of matter at SPS; a saturation of that 
quantity at the RHIC energies can be treated as an existence of the limiting Hagedorn temperature 
of hadronic matter, or maximal temperature of deconfinement T c . It is noteworthy that observation 
as for compatible values of the momentum-dependent APSD (f(p)) for collisions of different nuclei 
at top SPS and top AGS energies has been interpreted in Ref. jlUj as the universal properties 
of kinetic freeze-out in heavy ion collisions. Our study shows, however, that the pion APSD 
is approximately conserved value and so has no direct link to the freeze-out criterion and final 
thermodynamic parameters, being connected rather to the initial phase-space density of hadronic 
matter in A+A collisions [5]. 

A behavior of the entropy of thermal pions and measured pion multiplicities in central ra- 
pidity region vs energy demonstrates an anomalously high slope of an increase of the pion en- 
tropy/multiplicities at SPS energies compared to what takes place at the AGS and RHIC energies. 
This additional growth could be, probably, a manifestation of the QCD critical end point. The 
observed phenomenon can be caused by the dissipative effects that usually accompany phase transi- 
tions, such as an increase of the bulk viscosity, and also by peculiarities of pionic decays of a mesons 
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and other resonances with masses that are reduced, as compare to its vacuum values, in vicinity 
of the QCD CEP. At the RHIC energies there is no anomalous rise of pion entropy /multiplicities, 
apparently, because the crossover transition takes place far from the CEP and no additional degrees 
of freedom appear at that scale of energies: quarks and gluons were liberated at previous energy 
scale. 

We also analyze a behavior of the interferometry radii with energy in a context of the HBT 
puzzle. We show that if cubic power of the effective temperature of pion spectra grows with energy 
similar to the rapidity density then the interferometry volume is inversely proportional to the pion 
APSD. The behavior of the latter with collision energy is nearly constant starting from the high 
SPS energies because of the limiting Hagedorn temperature, T c , for hadronic matter. Roughly, 
the similarity between rise of effective temperature cubed and pion rapidity density takes place 
within the wide interval: from the lowest SPS energy to the highest RHIC energy. Therefore 
interferometry volume in Pb+Pb and Au+Au cental collisions is nearly constant, more precisely, it 
grows much slower (mostly due to an increase of longitudinal radius associated with total lifetime 
of the system) than rapidity density. At the same time, at each fixed energy the pion APSD, 
rapidity density and effective temperature of pion spectra are approximately conserved during the 
evolution [§], therefore the interferometry volume, that is the function of above values, is nearly 
the same as at the initial moment of hadronic evolution, if it were measured. It explains the 
experimental observations |231 135j that the interferometry volumes are changed only a little with 
energy for central collisions of the same nuclei and, at the same energy, they are proportional to 
the initial system extension in non-central collisions and central collisions of nuclei with different 
atomic numbers. The further experimental analysis of a such type of the correlations between the 
interferometry volumes, initial system sizes, multiplicities and slopes of the transverse spectra are 
still needed to clarify the picture. 

Summarizing this work, we point out that the available interferometry data not only do not con- 
tradict possible dramatic transformations of the matter in A+A collisions, as it is usually concluded 
[HHIISI], but being analyzed properly give deep insight into the physics of the phase transitions. The 
precise localization of the collision energy region where the transition to QGP at finite net baryonic 
densities happens is very important for understanding of the physics of deconfinement transition. 
At the moment no single experiment has collected data that by themselves show non-monotonic 
behavior of physical observables as a function of collision energy in the AGS-SPS energy domain. 
This can throw doubt upon an experimental significance of the observed APSD behavior which 
leads to the treatment of that as the result of a phase transition and critical end point. In view of 
this the Compressed Baryonic Matter (CBM) project 55 at the future Facility for Antiproton and 
Ion Research (FAIR) in Darmstadt is particularly important because it makes possible the system- 
atic studies of heavy ion collisions at beam energy range between 10 and 40 AGeV. Our results as 
for the localization of the QCD CEP indicate the importance of the experiments with relatively 
light nuclei at the top SPS energy and future RHIC energy scans in which the corresponding part 
of the (T — hb) plane can be reached. Also, a testing of our prediction of the APSD suturation 
at the top RHIC and LHC energies because of the limiting Hagedorn temperature, associated with 
the deconfinement temperature at zero net baryonic densities, is of a great importance. 
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Figure 1: The average phase-space density of all negative pions at mid-rapidity, (27r) 3 (f(y)), 
(circles, squares, stars and triangles) as function of cm. energy per nucleon in heavy ion central 
collisions. 
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Figure 2: The average phase-space density of all negative pions at mid-rapidity, (2tt) 3 (/(y)), as 
function of cm. energy per nucleon at the AGS-SPS energy domain. The error bars take into 
account uncertainties of the calculated APSD, see the text for details. 
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Figure 3: The average phase-space density of thermal ("direct") negative pions, (2tt) 3 (f(y)) th 
(rhombus), and the average phase-space density of negative pions at the stage of chemical freeze- 
out, (27r) 3 {f(y)) ch (crosses), as functions of cm. energy per nucleon in heavy ion central collisions. 
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Figure 4: The chemical potential of thermal ("direct") negative pions, /x t ^, (rhombus) as function 
of cm. energy per nucleon in heavy ion central collisions. 
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Figure 5: The interferometry volumes, Vi n t = (2tt) 3 / 2 RoRsRii (circles, squares, and stars) of 
negative pions at px — 0.06 — 0.07 GeV vs rapidity densities of the negative pions, dN n /dy, at 
mid-rapidity in heavy ion central collisions at different energies. 
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Figure 6: The ratio of rapidity densities of all negative pions to the corresponding interferometry 
volumes, (dN n /dy)/Vi n t, (circles, squares and stars) and the ratio of rapidity densities of negative 
thermal pions to their effective volumes, that is thermal densities nth, (crosses) vs cm. energies 
per nucleon in heavy ion central collisions. 
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Figure 7: The rapidity density of entropy for negative thermal pions, dS£ h /dy, (squares and stars) 
as function of cm. energy per nucleon in heavy ion central collisions. 
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Figure 8: The rapidity density of negative pions, dN n /dy, as function of cm. energy per nucleon 
in heavy ion central collisions. 
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